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Abstract 

We confront lattice QCD results on the transition from the hadronic phase to 
the quark-gluon plasma with hadron resonance gas and percolation models. We 
argue that for T <T C the equation of state derived from Monte-Carlo simulations 
of (2+1) quark-flavor QCD can be well described by a hadron resonance gas. We 
examine the quark mass dependence of the hadron spectrum on the lattice and 
discuss its description in terms of the MIT bag model. This is used to formulate a 
resonance gas model for arbitrary quark masses which can be compared to lattice 
calculations. We finally apply this model to analyze the quark mass dependence of 
the critical temperature obtained in lattice calculations. We show that the value 
of T c for different quark masses agrees with lines of constant energy density in a 
hadron resonance gas. For large quark masses a corresponding contribution from 
a glueball resonance gas is required. 



* Dedicated to Rolf Hagedorn 
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1 Introduction 



Long before lattice calculations provided first evidence [I] for critical behavior 
in strongly interacting matter it has been noticed [2] by Hagedorn that ordi- 
nary hadronic matter cannot persist as a hadronic resonance gas at arbitrary 
high temperatures and densities. This lead to the concept of the Hagedorn 
limiting temperature. With the formulation of QCD it has been suggested 
[3] that a phase transition to a new form of matter, the quark-gluon plasma, 
will occur. 

Two basic properties of hadrons were essential for developing the concept 
of a natural end for the era of ordinary hot and dense hadronic matter. In 
high energy experiments it had been observed that strongly interacting par- 
ticles produce a large number of new resonances. Moreover, hadrons have 
been known to be extended particles with a typical size of about 1 fm. As 
the average energy per particle increases at high temperatures copious par- 
ticle production will take place in a hadron gas and a dense equilibrated 
system will result from this. At high temperature, extended hadrons thus 
would start to "overlap" . This led to the expectation that some form of new 
physics has to occur under such conditions. The expected critical behavior 
has been analyzed in terms of various phenomenological models which incor- 
porate these basic features (resonance production =>■ Hagedorn's bootstrap 
model j2]; extended hadrons =>- percolation models In fact, many of the 
basic properties of the dense matter created today in heavy ion experiments 
can be understood quite well in terms of the thermodynamics of a hadronic 
resonance gas E] • 

With the formulation of Quantum Chromodynamics (QCD) as a theoret- 
ical framework for the strong interaction force among elementary particles it 
became clear that this "new physics" indeed meant a phase transition to a 
new phase of strongly interacting matter - the quark-gluon plasma (QGP) 
[3]. As QCD is an asymptotically free theory, the interaction vanishes loga- 
rithmically with increasing temperature, it has been expected that at least 
at very high temperatures the QGP would effectively behave like an ideal 
gas of quarks and gluons. Today we have a lot of information from numeri- 
cal calculations within the framework of lattice regularized QCD about the 
thermodynamics of hot and dense matter which give support to these ex- 
pectations. We know about the transition temperature to the QGP and the 
temperature dependence of basic bulk thermodynamic observables such as 
the energy density and the pressure [7j. In the coming years the increase in 
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numerical accuracy certainly will lead to modifications of the quantitative 
details of these results. However, already today they are sufficiently accurate 
to be confronted with theoretical and phenomenological models that provide 
a description of thermodynamics of strongly interacting matter. Recently, 
progress has been made to develop and link an improved perturbation the- 
ory of QCD with lattice data on the equation of state in the deconfined phase 
[S]. In this paper we analyze in how far the critical behavior can be under- 
stood in terms of the physical degrees of freedom of the confined phase, li.e. 
those of a hadronic resonance gas, and the intuitive percolation picture [5] . 

Quite distinct from the phenomenological approaches to the QCD phase 
transition are attempts to understand the thermodynamics of strongly inter- 
acting matter in terms of low energy effective theories, i.e. chiral perturba- 
tion theory JU| and effective chiral models [TT1 IT2] . The strength of these 
approaches is that they incorporate the correct symmetries of the QCD La- 
grangian and thus have a chance to predict the universal properties, e.g. the 
order of the phase transition, in the chiral limit of QCD. They, however, 
generally ignore the contributions of heavier resonances to the QCD thermo- 
dynamics which might be crucial for the transition to the plasma phase at 
non-vanishing values of the quark masses. 

Lattice calculations provide detailed information on the quark mass de- 
pendence of the transition to the QGP as well as to the hadron spectrum 
at zero temperature. In particular, we know that the transition temperature 
drops substantially when decreasing the quark mass from infinity (pure SU(3) 
gauge theory) to values close to the physical quark mass. This drop in the 
critical temperature can be understood at least qualitatively in terms of the 
relevant degrees of freedom in the low temperature phase. In the pure gauge 
limit this phase consists of rather heavy glueballs (mc>1.5 GeV [T31 EJ E3)- 
Quite a large temperature thus is needed to build up a sufficiently large den- 
sity of glueballs, which could lead to critical behavior (percolation [9 ). In the 
chiral limit, on the other hand, the low critical temperature can be addressed 
to the presence of light Goldstone particles, the pions, which can build up a 
large (energy) density already at rather low temperatures. Along with this 
decrease of the critical temperature goes an increase in the critical energy 
density expressed in units of the critical temperature, e c /T c 4 , by about an 
order of magnitude. This reflects the importance of new degrees of freedom 
in the presence of light quarks. However, at the same time the critical en- 
ergy density in physical units (GeV/fm 3 ) turns out to be almost quark mass 
independent. 
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In this paper we want to focus on these results. We will discuss in how 
far the quark mass dependence of the transition temperature found in lattice 
calculations is consistent with phenomenological models and what this tells 
us about the influence of the chiral sector of QCD on the transition tem- 
perature. In Section 2 we will briefly summarize the formulation of hadron 
thermodynamics in terms of a hadronic resonance gas. In section 3 we discuss 
the quark mass dependence of the hadron spectrum and give a phenomeno- 
logical parametrization motivated by the bag model. Predictions of these 
phenomenological approaches for the equation of state and the quark mass 
dependence of the transition temperature are then compared with lattice 
results in Section 4. Finally we give our conclusions in Section 5. 

2 Hadron resonance gas and the equation of 
state on the lattice 

The basic quantity required to verify thermodynamic properties of QCD is 
the partition function 1 Z(T,V). The grand canonical partition function is 
obtained as 



where H is the Hamiltonian of the system and (3 = 1/T is the inverse tem- 
perature. The confined phase of QCD we model as a non-interacting gas of 
resonances - the hadron resonance gas model. To do so we use as Hamilto- 
nian the sum of kinetic energies of relativistic Fermi and Bose particles of 
mass rrii. The main motivation of using this Hamiltonian is that it contains 
all relevant degrees of freedom of the confined, strongly interacting matter 
and implicitly includes interactions that result in resonance formation In 
addition this model was shown to provide a quite satisfactory description of 
particle production in heavy ion collisions [o^ lol llbj. 

With the above assumption on the dynamics the partition function can 
be calculated exactly and expressed as a sum over one-particle partition 
functions Z\ of all hadrons and resonances, 



1 We restrict our discussion to the case of vanishing chemical potential (vanishing net 
baryon number) and charge neutral systems. 



Z(T, V)=Tr[e-P B ] 



(1) 



\nZ(T,V) = J2^Zl(T,V). 



(2) 
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For particles of mass m ; and spin-isospin degeneracy factor the one- 
particle partition function Z\ is given by, 



\nZ}(T,V) = J£ J^dpp'rjHl + rje-^), (3) 

where = yjp 2 + m? is the particle energy and r] = — 1 for bosons and +1 
for fermions. 

Due to the factorization of the partition function in Eq. El the energy 
density and the pressure of the hadron resonance gas, 

e = £4 . P-Y.Pl (4) 

i i 

are also expressed as sums over single particle contributions e\ and P/, re- 
spectively. These are given by 



4 _ 9i v 2 * i ^k+i (PmY 



kpm 



(5) 



T 4 2vr 2 f-J 1 " k 

k=l 

a! s el ^^=ir* t(-vr i{ -^UkPm) (6) 

where i^i and K 2 are modified Bessel functions. 

Summing up in Eq. (0J) the contributions from experimentally known 
hadronic states, constitutes the resonance gas model for the thermodynam- 
ics of the low temperature phase of QCD. We take into account all mesonic 
and baryonic resonances with masses up to 1.8 GeV and 2.0 GeV, respec- 
tively. This amounts to 1026 resonances. The energy density obtained in this 
way starts rising rapidly at a temperature of about 160 MeV. In Fig. ^ we 
show the temperature dependence of the energy density e and the interaction 
measure A for the hadron resonance gas obtained from Eqs. (J1J) to (fB|). The 
model predictions are compared with Monte-Carlo results obtained ^7] on 
the lattice in (2+1) flavor QCD. Although it should be noted that the lattice 
calculations have not yet been performed with the correct quark mass spec- 
trum realized in nature the resonance gas model and the lattice data agree 
quite well. This indicates that for T < T c hadronic resonances are indeed the 
most important degrees of freedom in the confined phase. The energy density 
in the resonance gas reaches a value of 0.3 GeV/fm 3 at T ~ 155 MeV and 
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Figure 1: The left-hand figure shows the energy density e in units of T 4 
calculated on the lattice with (2+1) quark flavors as a function of the T/T c 
ratio. The vertical lines indicate the position of the critical temperature. 
The right-hand figure represents the corresponding results for the interaction 
measure (e — 3P)/T 4 . The full-lines are the results of the hadron resonance 
gas model that accounts for all mesonic and baryonic resonances. 

1 GeV/fm 3 already at T ~ 180 MeV. This is in good agreement with lattice 
calculations, which find a critical energy density of about 0.7 GeV/fm 3 at 
T c ~ 170 MeV ^7j. For comparison we note that a simple pion gas would 
only lead to an energy density of about 0.1 GeV/fm 3 at this temperature. 
This suggests that a more quantitative comparison between numerical results 
obtained from lattice calculations and the resonance gas model might indeed 
be meaningful. 



3 Hadron spectrum in heavy quark— mass limit 

In order to use the resonance gas model for further comparison with lattice 
results we should take into account that lattice calculations are generally 
performed with quark masses heavier than those realized in nature. In fact, 
we should take advantage of this by comparing lattice results obtained for 
different quark masses with resonance gas model calculations based on a 
modified, quark mass dependent, resonance spectrum. 

Rather than converting the bare quark masses used in lattice calculation 
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into a renormalized mass it is much more convenient to use directly the pion 
mass (m 7T ~ y/m q ), i.e. the mass of the Goldstone particle, as a control 
parameter for the quark mass dependence of the hadron spectrum. For our 
thermodynamic considerations we need, at present, not be concerned with 
the detailed structure of the hadron spectrum in the light quark mass chiral 
limit. We rather want to extract information on the gross features of the 
quark mass dependence of a large set of resonances. In order to study the 
quark mass dependence of hadron masses in the intermediate region between 
the chiral and heavy quark mass limits we adopt here an approach that is 
based on the Hamiltonian of the MIT bag-model ^H]- Although, in the 
original formulation this Hamiltonian breaks explicitly chiral symmetry and 
implies non-conservation of the axial-vector current it still provides a sat- 
isfactory description of the hadron mass spectrum that can be used for our 
thermodynamic considerations. 

In the limit of a static, spherical cavity the energy of the bag of radius R 
is given by 

E = E V + E + E K + E M + E E . (7) 

The first two terms are due to quantum fluctuations and are assumed to 
depend only on the bag radius. The volume and the zero-point energy terms 
have a generic form 

E v = ~nBR 3 , £o = -§ , (8) 

where B is the bag constant and Zq is a phenomenological parameter at- 
tributed to the surface energy. 

The quarks inside the bag contribute with their kinetic and rest energy. 
Assuming N quarks of mass rrii the quark kinetic energy is determined from 

1 N 

e k = + i^R) 2 ] 1/2 , (9) 
H i=l 

where x^im^ R) enters the expression on the frequency lo = [x 2 + (mR) 2 ] 1 / 2 / R 
of the lowest quark mode and is obtained [TSj as the smallest positive root 
of the following equation 
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tan(xj) = j . (10) 

1 — rriiR — Jxf + (rriiR) 2 

The last two terms in Eq. (JJJ) represent the color-magnetic and electric 
interaction of quarks. It is described by the exchange of a single gluon be- 
tween two quarks inside the bag. The color electric energy was found in [IB] 
to be numerical small and will be neglected in our further discussion. The 
color magnetic exchange term is given by 

c , ^ M {rriiR, rrijR) 

E M = Ska c }_^ — J - — (a ■ <Tj) . (11) 

Here a c is the strong coupling constant and k — 1 for baryons and 2 for 
mesons. For a given spin configuration of the bag the scalar spin product 
in Eq. (JTTJ) can easily be calculated. The function M(x, y) depends on the 
quark modes magnetic moment and is described in detail in ^HJ- For small 
x < 1 it shows a linear dependence on the argument with M(0, 0) = 0.175. 

The dependence of the energy on the bag radius can be eliminated by 
the condition that the quark and gluon field pressure balance the external 
vacuum pressure. For a static spherical bag this condition is equivalent to 
minimizing E with respect to R. The true radius R of the bag thus is 
determined from the condition BE / OR = and the hadron bag mass is then 
obtained from Eq. (j7J) with R = R Q . 

To extract the physical mass spectrum from the MIT bag model one 
still needs to fix the set of five parameters that determine the bag energy. 
Following the original fit to experimental data made in ^HJ we take: B 1 ^ = 
0.145 GeV, Z Q = 1.84 and a c = 0.55. These parameters together with 
Tn u = ma = and m s = 0.279 GeV provide a quite satisfactory description 
of hadron masses belonging to the octet and decuplet of baryons and the 
octet of vector mesons. Of course, the model fails to describe the details of 
the chiral limit and, in particular, it leads to a too large value of the pion 
mass that with the above set of parameters is m n = 0.28 GeV. Nonetheless, 
the accuracy of the bag model will be sufficient for our purpose. 

The MIT bag model provides an explicit dependence of hadron masses 
on the constituent quark mass. This dependence is entirely determined by 
the kinetic and magnetic energy of the quarks. To compare bag model calcu- 
lations with lattice calculations, which do not provide values for constituent 
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Figure 2: Dependence of different hadron masses rrih on the pion mass m^. 
Both and are expressed in the units of the string tension y/a. Curves 
are the MIT bag model results (see text for details). The filled circles rep- 
resent the PC-PACS lattice results from j2Bj- The filled diamonds are the 
Nf = 3 whereas the open-diamonds are Nf = 2 flavor results from |2"2"] . 
The filled-boxes are quenched QCD results [21] • All other points are from 
reference . Both the lattice data and the bag model results are shifted in 
m^-direction by a constant factors indicated in the figure. 



quark masses, it is best to express the quark mass dependence in terms of 
the pion mass, which is most sensitive to changes of the quark masses. In 
Fig. |21 we show the resulting dependence of different hadron masses on the 
pion mass with the bag parameters described above but with varying m u . 
The masses are expressed in units of the square root of the string tension 
for which we use yfa = 420 MeV. The model predictions are compared with 
recent lattice data on hadron masses calculated for different current quark 
masses |221 ESI EH ES! The MIT bag model is seen in Fig. to describe lat- 
tice results quite well. This is particularly the case for larger quark masses 
such that > y/a. For < yfa the deviations of the model from the 
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lattice results are quite apparent. As mentioned this is, of course, mainly 
due to the well known limitations of the bag model when approaching the 
chiral limit. 

For large quark masses the bag model description of hadron masses re- 
produces the naive parton model picture and consequently all hadron masses 
are almost linearly increasing with the pion mass as seen in Fig. |2] . This is 
to be expected as in this case the energy of the bag is entirely determined by 
the quark rest mass. As seen in Fig. El the slope increases with the number 
of non-strange constituent quarks inside the bag. Consequently, the slops of 
(H*,H) and (K,K*) or (£*,A) and p coincide at large m n . 

In order to formulate a resonance gas model for arbitrary quark masses we 
need to know the quark mass dependence of much more resonances than the 
few hadronic states shown in Fig. |21 We thus looked for a phenomenological 
parametrization of the quark mass dependence of resonances, expressed in 
terms of the pion mass. Fig. |21 suggests that already at intermediate values 
of the quark mass, rn^y^a, this dependence is dominated by the quark 
rest mass and does not depend much on the hadronic quantum numbers. 
This suggests that a common parametrization of all hadronic states, which 
is consistent with the naive parton model picture for large quark masses and 
reproduces the experimental values of hadronic states in the light quark mass 
limit is sufficient for our thermodynamic considerations. To incorporate these 
features we use the ansatz, 



which provides a good description of the MIT bag model result for non- 
strange hadron masses calculated for different values of m n . Here x = 
m n/\/&, n^o = Trihadron/y/o', N u is the number of light quarks inside the 
hadron (N u = 2 for mesons N u = 3 for baryons) and a = (0.42 GeV) 2 is the 
string tension. 

The parameters appearing in Eq. (|12p were optimized such that they 
reproduce the MIT bag model results for the m^-dependence of the p vector 
meson mass and are summarized in Table 1. In the mass regime shown in 
Fig. |21 Eq. (|T2~j) reproduces the quark mass dependence of all non-strange 
hadron masses obtained from the bag model within a relative error of <6%. 

In the following we will use Eq. (|T2"j) to formulate a hadron resonance gas 
model with varying quark masses. It will then be compared with lattice cal- 



M(x) 



~ N u aix + 



1 + a 2 x + a 3 x 2 + a 4 x 3 + a 5 x 4 



(12) 
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CLi 


a 2 


a 3 


04 


a 5 


0.51± 0.1 


aiN u 
m 


0.115± 0.02 


-0.0223± 0.008 


0.0028± 0.0015 



Table 1: Parameters entering the interpolation formula for non-strange 
hadron masses given in Eq. (12). 



culation of QCD thermodynamics. We will test, in particular, in how far this 
model can provide a quantitative description of the transition temperature 
obtained on the lattice for different quark masses. 

4 Quark mass dependence of the QCD tran- 
sition 

We want to confront here the resonance gas model developed in the previous 
section with lattice results on the quark mass dependence of the QCD tran- 
sition temperature and use it to learn about the critical conditions near de- 
confinement. Lattice calculations suggest that this transition is a true phase 
transition only in small quark mass intervals in the light and heavy quark 
mass regime, respectively. In a broad intermediate regime, in which the pion 
mass changes by more than an order of magnitude, the transition is not re- 
lated to any singular behavior of the QCD partition function. Nonetheless, it 
still is well localized and is characterized by rapid changes of thermodynamic 
quantities in a narrow temperature interval. The transition temperature thus 
is well defined and is determined in lattice calculations through the location 
of maxima in response functions such as the chiral susceptibility. A collec- 
tion of transition temperatures obtained in calculations with 2 and 3 quark 
flavors with degenerate masses is shown in Fig. El The main feature of the 
numerical results which we want to explore here is that the transition tem- 
perature varies rather slowly with the quark mass. In Ref. J7j the almost 
linear behavior has been described by the fit, 
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Figure 3: The transition temperature in 2 (filled squares) and 3 (circles) 
flavor QCD versus mps/y/o 1 using an improved staggered fermion action (p4- 
action). Also shown are results for 2-flavor QCD obtained with the standard 
staggered fermion action (open squares). The dashed band indicates the 
uncertainty on T c j ^fo in the quenched limit. The straight line is the fit 
given in Eq. EH 

which also is shown in Fig. El For pion masses rrtps ~ (6 — 7)y/a ~ 2.5 GeV 
the transition temperature reaches the pure gauge value, T c /*Ja ~ 0.632(2) 

E3- 

We note that all numerical results shown in Fig. El do correspond to quark 
mass values in the crossover regime. Also the resonance gas model formulated 
in the previous section does not lead to a true phase transition. We thus may 
ask what the conditions in a hadron gas are that trigger the transition to the 
plasma phase. Using the hadron gas with a quark mass dependent hadron 
mass spectrum and including the same set of 1026 resonances which have been 
included in other phenomenological calculations |H] we have constructed 
resonance gas models for 2 and 3 flavor QCD, respectively. In the former case 
we eliminate all states containing strange quarks whereas in the latter case 
we assigned to meson states containing strange particles the corresponding 
masses of non-strange particles, e.g. kaons have been replaced by pions etc. 
With these resonance gas models we have calculated the energy density at 
the transition temperature. We use T c = 175 (15) MeV for 2-flavor QCD and 
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T c = 155 (15) MeV for 3-flavor QCD, respectively. For the energy densities 
at the transition point we then find 



This is in good agreement with the lattice result, e/T c 4 = (6 ± 2) quoted in 
|2H1 as an average for the 2 and 3-flavor energy densities. In fact, as can 
be seen from Fig. 5 in Ref. [2H] the difference in e/T c 4 in the lattice results is 
of similar magnitude as we found here from the resonance gas model. The 
lattice results for 2 and 3-flavor QCD thus suggest that the conditions at the 
transition point are well described by a resonance gas. For comparison we 
also note that in the 2-flavor case a pion gas does contribute only about 20% 
to this energy density 2 and also a gas build up from the 20 lowest resonances 
would give rise only to about half the critical energy density, i.e. e/T 4 ~ 1.9. 

Although the lattice results allow, at present, only to determine the criti- 
cal energy density within a factor (2-3) it is striking that the transition occurs 
at similar values of the energy density in QCD with light quarks as well as in 
the pure gauge theory, although the transition temperature shifts by about 
40% and e/T c 4 differs by an order of magnitude. It thus has been suggested 
that for arbitrary quark masses the transition occurs at roughly constant en- 
ergy density. Such an assumption is, in fact, supported by our resonance gas 
model constructed in the previous section for arbitrary values of the quark 
masses. In Fig. |1] we show lines of constant energy density calculated in the 
resonance gas model and compare these to the transition temperatures ob- 
tained in lattice calculations. As can be seen the agreement is quite good up 
to masses, mps — 3 \fo or mps — 1.2 GeV. The reason for the deviations 
at larger values of the quark mass, of course, is due to the fact that we have 
neglected so far completely the glueball sector in our considerations. 

When the lightest hadron mass becomes comparable to typical glueball 
masses, also the glueball sector will start to contribute a significant fraction 
to the energy density. Using the set of 15 different glueball states so far 
identified in lattice calculations ^3] we have calculated their contribution to 
the total energy density. At vripsj \fo — 6.5 they contribute as much as the 
entire hadronic sector. However, as can be seen in Fig. El the contribution 
of these 15 states only leads to a small shift in the lines of constant energy 
density. Similar to the hadronic resonance gas for small quark masses where 

2 For massless pions we have e/T 4 = (n? — l)7r 2 /30 ~ 1. 




(14) 
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Figure 4: The transition temperature vs. pion mass obtained in lattice calcu- 
lations and lines of constant energy density calculated in a resonance gas model. 
The left hand figure shows a comparison of constant energy density lines at 1.2 
(upper), 0.8(middle) and 0.4(lower) GeV/fm 3 with lattice results for 2-flavor QCD 
obtained with improved staggered %Tj\ as well as improved Wilson [201 l2T] 
fermion formulations. T c as well as mpg are expressed in terms of the correspond- 
ing vector meson mass. The right hand figure shows results for 2 and 3 flavor QCD 
compared to lines of constant energy density of 0.8 GeV/fm 3 . Here T c and mps 
are expressed in units of y/a. For a detailed description see text. 



the 20 low-lying states only contribute 50% of the total energy density one 
has to expect that also in the large quark mass limit further glueball states, 
which have so far not been identified, will contribute to the thermodynamics. 
Further support for this comes from a calculation of the energy density of 
the 15 known glueball states at the transition temperature of the pure gauge 
theory, T = 0.63^. For this we obtain e(T = 0.63^) ~= .06 GeV/fm 3 
or equivalently e/T c 4 ~ 0.1, which is about 20% of the overall energy density 
at T c . The contribution of the 15 glueball states thus does not seem to 
be sufficient. In fact, the transition temperature in d- dimensional SU(N C ) 
gauge theories is well understood in terms of the critical temperature of string 
models, T c j y/a = \j3/ir(d — 2), which also is due to an exponentially rising 
"mass" spectrum for string excitations [2*Uj . 

It is conceivable that extending the glueball mass spectrum to all higher 
excited states will improve the results shown in Fig. |S| On the other hand 
one also should stress that the glueball states used in our calculations were 
obtained in quenched QCD and at zero temperature. There are indications 
from lattice calculations that glueball masses could be modified substantially 
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Figure 5: The transition temperature in 3-flavor QCD compared to lines 
of constant energy density (e = 0.8 GeV/fm 3 ) in a hadronic resonance gas 
(upper curve), a hadronic resonance gas with 15 glueball states added (middle 
curve) and a hadronic resonance gas with 15 glueball states with a 40% 
reduced mass (lower curve). 



in the presence of dynamical quarks jT3] as well as at finite temperature ^3] . 
The analysis of glueball states at high temperature [Hj suggests that their 
masses can drop by ~ (20 — 40)%. As all glueballs are heavy on the temper- 
ature scale of interest, shifts in their masses influence the thermodynamics 
much more strongly than in the light quark mass regime where the lowest 
state has already a mass which is of the order of the transition temperature. 
In fact, we find that taking into account a possible decrease of the glueball 
masses close to T c seems to be more important than adding further heavy 
states to the spectrum. We thus have included a possible reduction of glue- 
ball masses in the equation of state. The resulting T c with this modification 
is also shown in Fig. |5J Decreasing the glueball masses, increases the ther- 
mal phase space available for particles, thus consequently the temperature 
required to get e = 0.8 GeV/fm 3 is decreasing. As can be seen in Fig. El a 
reduction of glueball masses by 40% is sufficient to reproduce lattice results 
in the whole m n range. However, to make this comparison more precise it 
clearly is important to get a more detailed understanding of the glueball 
sector in the future. 
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5 Conclusions 



In this paper we have analyzed lattice results on QCD thermodynamics using 
a phenomenological hadron resonance gas model. We have shown that close 
to the chiral limit and for T <T C the equation of state derived on the lattice 
is quantitatively well described by the resonance gas. 

The hadron resonance gas partition function is also shown to be suitable 
to describe lattice results for finite quark masses and varying number of 
flavors. One needs, however, to implement the quark mass dependence of 
the hadron spectrum and for large values of the quark mass the glueball 
degrees of freedom have to be taken into account as they start playing an 
important role. We have shown that, away from the chiral limit region, 
the quark mass dependence arising from MIT bag model agrees quite well 
with the hadron mass spectrum calculated on the lattice. We find that the 
transition temperatures obtained in lattice calculations at different values of 
the quark mass are well described by lines of constant energy density in a 
resonance gas model. For moderate values of the quark masses the predictions 
of the hadron resonance model coincide with lattice calculations. However, 
for heavy quark masses this agreement could be only achieved by including 
additional heavy glueball states or allowing for a reduction of glueball masses 
close to the transition temperature by about 40%. 

Our results can be considered as an indication that thermodynamics in 
the vicinity of deconfinement is indeed driven by the higher excited hadronic 
states. This finding can give additional support for previous phenomenolog- 
ical applications of the resonance gas partition function in the description 
of particle production in heavy ion collisions. Our discussion of the critical 
temperature and its quark mass dependence also indicates that deconfine- 
ment in QCD to large extend is density driven. It would be interesting to 
see to what extent the lines of constant energy density of the generalized 
hadron resonance gas can be related to correspondingly generalized percola- 
tion models. 
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